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ISOTHERAAAL ELASTOHYDRODYNAMIC LUBRICATION OF POINT CONTACTS 

III - FULLY FLOODED RESULTS* 
by Bernard J. Hamrock and Duncan Dowson* 

Lewis Research Center 

SUMMARY 

The theory developed by the authors in an earlier publication was used to investigate 
the influence of the ellipticity parameter and the dimensionless speed U, load W, and 
material parameters on minimum film thickness. The ellipticity parameter was varied 
from 1 (a ball-on-plate configuration) to 8 (a configuration approaching a line contact). 
The dimensionless speed parameter was varied over a range of nearly two orders of 
magnitude. And the dimensionless load parameter was varied over a range of one order 
of magnitude. Conditions corresponding to the use of solid materials of bronze, steel, 
and silicon nitride and lubricants of paraffinic and naphthenic mineral oils were consid- 
ered in obtaining the exponent on the dimensionless material parameter. Thirty -four 
different cases were used in obtaining the minimum -film -thickness formula 

®min = 3,63 U0.68 G 0.49 W -0.07 3 (1 _ e ~0. 68kj 
A simplified expression for the ellipticity parameter k was found, where 

/R \ 0 * 64 
k=1.03(-X\ 

VV 


R /R being the radius of curvature ratio, 

y ^ 

Contour plots are also shown that indicate in detail the pressure spike and the two 
side lobes in which the minimum film thickness occurs. These theoretical solutions of 
film thickness have all the essential features of previously reported experimental ob- 
servations based upon optical interferometry. 

* Presented at Joint Lubrication Conference cosponsored by the American Society of 
Mechanical Engineers and the American Society of Lubrication Engineers, Boston, Mas- 
sachusetts, October 5-7, 1976. 

^ Professor of Mechanical Engineering, University of Leeds, Leeds, England. 



INTRODUCTION 


Only in recent years has the complete theoretical solution of the isothermal elasto- 
hydrodynamic lubrication (EHL) of point contacts been successfully analyzed. The anal- 
ysis requires the simultaneous solution of the elasticity and Reynolds equations . The 
authors’ approach to the theoretical solution has been presented in two previous publica- 
tions (refs. 1 and 2) . The first of these publications (ref. 1) presents the elasticity 
model in which the conjunction is divided into equal rectangular areas with a uniform 
pressure applied over each area. The second (ref. 2) gives the complete approach to 
solving the elastohydrodynamic lubrication problem for point contacts. 

The most important practical aspect of the EHL point-contact theory (ref. 2) is the 
determination of the minimum film thickness within the contact. That is, maintaining a 
fluid film thickness of adequate magnitude is extremely important to the operation of 
some machine elements. In the present report, only the results from the theory given 
in references 1 and 2 are presented. In the results the influence of contact geometry as 
expressed in the ellipticity parameter and the dimensionless speed, load, and material 
parameters on minimum film thickness is investigated for a conjunction fully immersed 
in lubricant (i.e. , fully flooded) . The ellipticity results have been presented in refer- 
ence 3 by the authors, but in the present report additional solutions are presented and 
the accuracy of the presentation is improved. The corrected ellipticity results are pre- 
sented in this report. In the numerical work the ellipticity parameter was varied from 
1 (a ball -on-plate configuration) to 8 (a configuration approaching a line contact) . The 
dimensionless speed and load parameters were varied over ranges of about two and one 
orders of magnitude, respectively. Conditions equivalent to using solid materials of 
bronze, steel, and silicon nitride and lubricants of paraffinic and naphthenic mineral 
oils were considered in obtaining the exponent on the dimensionless material parameter. 
Thirty -four different cases were used in obtaining the fully flooded film-thickness for- 
mulas. A fully flooded condition is said to exist when the inlet distance of the conjunc- 
tion ceases to influence, in any significant way, the minimum film thickness. The inlet 
distance of the conjunction is defined as the distance from the center of the contact to the 
edge of the computing area. Contour plots are also shown that indicate in detail the 
pressure spike and the two side lobes in which the minimum film thickness occurs. 


SYMBOLS 


a semimajor axis of contact ellipse 
b semiminor axis of contact ellipse 
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H min 


H min 
^min, L 


x, X,X 

y, Y,Y 



modulus of elasticity 
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normal applied load 

dimensionless material parameter, E'/p^y as 
dimensionless film thickness, h/R x 

dimensionless central film thickness obtained from EHL point -contact theory 
(ref. 2) 

dimensionless central film thickness obtained from least-square fit of data 

dimensionless minimum film thickness obtained from EHL point-contact 
theory (ref. 2) 

dimensionless minimum film thickness obtained from least-square fit of data 

dimensionless minimum film thickness for line contact 

film thickness 

ellipticity parameter, a/b 

dimensionless pressure, p/E’ 

pressure 

asymptotic isoviscous pressure 
effective radius 
radius of curvature 

dimensionless speed parameter, utj 0 /E'R x 
surface velocity in x-direction 

o 

dimensionless load parameter, F/E'R x 
coordinate systems defined in report 
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a pressure -viscosity coefficient 
7J o atmospheric viscosity 
v Poisson’s ratio 

Subscripts: 

A solid A 

B solid B 

x,y coordinates system defined in report 

DIMENSIONLESS GROUPING 

From the variables of the numerical analysis (ref. 2) the following dimensionless 
groups can be defined: 

(1) Dimensionless film thickness 


(2) Ellipticity parameter 




(3) Dimensionless speed parameter 


E'R X 

(4) Dimensionless load parameter 



(5) Dimensionless material parameter 

a— 21- 

Piv , as 


( 1 ) 


( 2 ) 


(3) 


(4) 


( 5 ) 
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where p. is the asymptotic isoviscous pressure obtained from Roelands (ref. 4). 

IV J MrS 

The asymptotic isoviscous pressure can be approximated by the inverse of the pressure 
viscosity coefficient (p. « 1/a). 

XV ^ d o 

The dimensionless film thickness can be written as 


H=f(k,U,W,G) (6) 

The most important practical aspect of the EHL point-contact theory developed in refer- 
ence 2 is the determination of the minimum film thickness within the conjunction. 
Therefore, in the fully flooded results presented herein the dimensionless parameters 
(k, U, W, and G) will be varied and the effect upon minimum film thickness will be 
studied . 


Effect of Ellipticity of Elastic Conjunction 


The ellipticity parameter k is a function of the radii of curvature of the solids only 
( r Ax> r Bx> r Ay ’ and r By^ * The radii oi curvature in the x-direction for both solids A 
and B are used in defining the dimensionless speed and load parameters. Therefore, 
only the radius of curvature of solid B in the y-direction was changed in varying the el- 
lipticity parameter from 1 (a ball on -plate configuration) to 8 (a configuration approach- 
ing a line contact). In doing this the dimensionless speed U, load W, and material G 
parameters were held constant at the following values: 


U = 0.1683X10' 11 ^ 
W = 0.1106X10" 6 
G = 4522 


(7) 


Care was taken to ensure that the highest ellipticity parameter (k = 8) case was in the 
elastic region as defined by Dowson-Higginson (p. 101, ref. 5). For ellipticity param- 
eters less than 8 the results move farther into the elastic region. 

Table I gives 10 values of the ellipticity parameter k and the corresponding mini- 
mum film thickness H min as obtained from the EHL point-contact theory (ref. 2) . 
Having these 10 pairs of data, the object was to determine an equation that describes how 
the ellipticity parameter affects the minimum film thickness. The general form of this 
equation can be written as 
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( 8 ) 



TJ 

min 


H 


min, Ly 


= A e 


Bk 


A least-square exponential curve fit to the 10 pairs of data points 


I 1 ~ 


H 


min 


H min,L/. 


where i = 1 , . . . , 10 


was used in obtaining values for A and B in equation (8) . Besides a least-square fit a 

o 2 

coefficient of determination r was obtained . The value of r reflects the fit of the 

data to the resulting equation: 1 being a perfect fit, and zero the worst possible fit. The 

minimum film thickness for a line contact H min L used in equation (8) was determined 

by finding the H . T that gives a coefficient of determination closest to 1. This value 
mm, 

of H . t turned out to be 7.082x10 with a corresponding coefficient of determination 
min, J-i 

of 0.9990, which is an excellent fit. Furthermore, the values of A and B inequa- 
tion (8) as obtained from the least-square fit are 


A = 0.9966 * 1.00 (9) 

B = -0.6752 « -0.68 (10) 


From equations (8), (9), and (10) the following equation can be written, which shows 
the effect of ellipticity parameter on minimum film thickness: 




(XI) 


It is most significant that A turned out to be 0.9966, or approximately 1.00, since as 
k — 0, H m ^ n —0. Therefore, even though the smallest value of k used in obtaining 
equation (11) was unity, it would seem that equation (11) could be applied to smaller 
values since in the limiting case (k —0) equation (11) satisfies the physical intuition. 

For the other extreme of large k, a line -contact situation is approached, and the agree- 
ment with existing results is again good. From Dowson and Higginson (ref. 5) the line- 
contact minimum thickness for the dimensionless parameters given in equation (7) is 
7.720X10" 6 . Compare this with 7.082x10"® from the present results. The difference of 
9 percent could well be the result of Dowson and Higginson (ref. 5) using an exponential 
pressure -viscosity relation instead of the Roelands (ref. 4) formulation used in the pres- 
ent work. As was pointed out in the closure of reference 2 (answering a query by 
P. M. Ku) the Roelands formula (ref. 4) suppresses the pressure spike somewhat and 
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also results in a smaller film thickness. 

Substituting equations (9) and (10) into equation (8) gives H min the dimensionless 
minim um film thickness obtained from the least -square formulation. The H min for 
10 values of the ellipticity parameter are given in table I. The percentage difference 
between the mi nim um film thickness obtained from EHL point -contact theory H m ^ n and 
the minimum film thickness obtained from the least-square fit equation H m - n is ex- 
pressed as 


D l = 


<H _ H 

min min i 


10C 


H 


min 


( 12 ) 


Note that in table I, the magnitude of is within ±3 percent for all cases. 

Figures 1(a) and (b) give contour plots of dimensionless pressure for two extreme 
values of the ellipticity parameter k, 8 and 1.25. In these and all contour plots to be 
presented the + symbol indicates the center of the Hertzian contact. Note that, because 
of the dimensionless representation of the X- and Y -coordinates, the actual Hertzian 
contact ellipse becomes a circle regardless of the value of the ellipticity parameter . 

The Hertzian contact circle is shown in each figure by asterisks. At the top of each 
figure the contour labels and the corresponding values of dimensionless pressure are 
given. The inlet region is to the left and the exit region is to the right. 

For an ellipticity parameter of 8 the maximum pressure is near the center of the 
contact; and, even though the conditions are in the elastic region, no pressure spike 
occurs. The pressure gradient at the exit end of the conjunction is much larger than 
that in the inlet region. For an ellipticity parameter of 1.25 a pressure spike is visible 
at the exit end of the contact. 

Figures 2(a) and (b) show two contour plots of film thickness when the ellipticity pa- 
rameter k is 8 and 1.25, respectively. For an ellipticity parameter of 8 the minimum 
film thickness is directly behind the axial center of the contact. For an ellipticity pa- 
rameter of 1.25, two minimum-film -thickness regions occur to the sides and nearer the 
edge of the Hertzian circle. These results - showing the two ’’side lobes” in which 
minimum-film-thickness areas occur - produce all the essential features of previously 
reported experimental observations based upon optical interferometry (ref. 6). 

Figures 3(a) and (b) show the variation of pressure and film thickness, respectively, 
in the X-direction close to the midplane of the contact for three values of the ellipticity 
parameter. As has been true for all the ellipticity parameter results presented, the 
values of the dimensionless speed, load, and material parameters were held fixed as per 
equation (7) . In figure 3(a) we find that for k = 6 no pressure spike occurs, but this 
well-known feature of theoretical solutions to the EHL problem is evident for k = 2. 5 
and k = 1.25. 


7 



In figure 3(b) for k = 1.25 the central region is not parallel with the X-axis. The 
reason is probably that compressibility effects are considered in the theory (ref. 2). 
That is, when compressibility is considered, the film thickness in the center is reduced 
by the amount that the fluid volume decreases at high pressure (ref. 5) . 


INFLUENCE OF SPEED 

By changing only the surface velocity in the x-direction u, the dimensionless speed 
parameter U (eq. (3)) changes, but the other dimensionless parameters (k, W, and G) 
remain constant. The values at which these dimensionless parameters were held con- 
stant in the calculations are 


W = 0.7371X10" 6 ^ 
G = 4522 
k = 6 > 


(13) 


Table II gives the dimensionless speed parameter U and the corresponding mini- 
mum film thickness H ^ as obtained from the EHL point-contact theory (ref. 2). 
There are 15 different values of the dimensionless speed parameter covering nearly two 
orders of magnitude. Having these 15 pairs of data, the objective is to determine an 
equation that describes how the dimensionless speed affects the minimum film thickness. 
The general form of this equation can be written as 


H min = 


(14) 


By applying a least-square power fit to the 15 pairs of data (Up H min p where 
i = 1, . . . , 15) the values of I and J were found to be 


I = 560.18 


(15) 


J = 0.67542 » 0.68 (16) 

o 

The coefficient of determination r for these results was excellent at 0.9998. Substitut- 
ing equations (15) and (16) into equation (14) gives the values of H min shown in table II. 
The percentage difference between the minimum film thickness obtained from the 
EHL point -contact theory H m ^ n and the minimum film thickness obtained from the 
least-square fit H min is expressed in equation (12) and given in table H. Note that the 
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variation of Dj is less than ±2 percent. 

From equations (14) and (16) the effect of dimensionless speed on dimensionless 
minimum film thickness can be written as 


H . cc u 
min 


0.68 


( 17 ) 


Figures 4(a) and (b) give contour plots of dimensionless pressure for two extreme 
values of the dimensionless speed parameter U, 0.8416x10"^ and 0. 5050x10"^. In 
the low-speed case (U = 0. 8416xlO - ^) there is a pressure spike at the exit end of the 
contact. In the high-speed case (U = 0. 5050xl0~ 10 ) no pressure spike occurs. Note 
from figures 4(a) and (b) that the pressure in the inlet region is higher at high speeds 
than at low speeds . 

Figures 5(a) and (b) show contour plots of dimensionless film thickness when the di- 
mensionless speed parameter is 0. 8416xl0~ 1 ^ and 0. 5050X10 -10 , respectively. In fig- 
ure 5(a) the minimum film thickness appears close to the Hertzian circle and off to the 
side. In figure 5(b) the minimum -film-thickness area appears between the center of the 
contact and the Hertzian circle . 

Figures 6(a) and (b) show the variation of pressure and film thickness, respectively, 
on the X-axis near the midplane of the conjunction for three values of dimensionless 
speed parameter. As has been true for all the speed results, the values of the dimen- 
sionless load, material, and ellipticity parameters are held fixed as described in equa- 
tion (13). In figure 6(a) the dashed line corresponds to the Hertzian pressure distribu- 
tion. It can be seen from figure 6(a) that the pressure in the inlet region is higher for 
the high speed (U = 0. 5050x10"^) profile. For U = 0.8416x10”^ the pressure spike 
originates very near the trailing edge of the Hertzian pressure distribution, and as the 
speed increases the pressure spike moves upstream. For line contacts, Dowson and 
Higginson (ref. 5) found results similar to those shown in figure 6(a). 

The typical elastohydrodynamic film shape with an essentially parallel section in 
the central region is shown in figure 6(b) . Also, there is a considerable change in film 
thickness as the dimensionless speed is changed, as indicated by equation (17). This il- 
lustrates most clearly, the dominant effect of the dimensionless speed parameter upon the 
minimum film thickness in elastohydrodynamic contacts. 


INFLUENCE OF LOAD 

Changing only the normal applied load F in equation (4) causes the dimensionless 
load parameter W to change while the remaining dimensionless parameters (k, U, and 
G) remain constant. The values at which these parameters were held constant are 
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U = 0. 1683x10 ~ 11 ''| 
G = 4522 


( 18 ) 


k = 6 


Table in gives the dimensionless load parameter and the corresponding minimum 
film thickness H min as obtained from the EHL point-contact theory (ref. 2). There 
are eight different values of the dimensionless load parameter covering more than an 
order of magnitude . Having these eight pairs of data, the objective was to determine an 
equation that describes how the dimensionless load affects the minimum film thickness. 
The general form of this equation can be written as 


H 


mm 


KW 


(19) 


By applying a least-square power fit to the eight pairs of data (W. , H . ., where 

JL I1JUL1I j X 

i = 1, . . ., 8) the values of K and L were found to be 

K = 2. 1592X10 “ 6 (20) 

L = -0.072924 » -0.073 (21) 

2 

The coefficient of determination r for these results was 0.9260, which was good but 
was the lowest obtained in deriving the minimum -film -thickness equation (eq. (28)). 
Substituting equations (20) and (21) into equation (19) gives the values of H m ^ n shown in 
table HI. The percentage difference Dj between the minimum film thickness obtained 
from the EHL point-contact theory H m ^ n and the minimum film thickness obtained from 
the least-square fit equation H min is expressed in equation (12) and given in table HI. 

In table IH the variation of is within ±3 percent in all cases. 

From equations (19) and (21) the effect of load on minimum film thickness can be 
written as 


H min 


oc w 


-0.073 


( 22 ) 


Figures 7(a) and (b) give contour plots of dimensionless pressure for the two extreme 
values of dimensionless load parameter W that were investigated, 0.1106X10 -6 and 
0.1290x10”^. No pressure spike is evident at the lowest load, but at the highest load a 
pressure spike is visible. 

Contour plots of dimensionless film thickness for the same two values of dimension- 
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II 


less load parameter as given in figure 7 are shown in figure 8. In figure 8(a) for the 
low-load case (W = 0. 1106x10"®) the minimum film thickness occurs directly behind the 
center of the contact. In figure 8(b) for the high -load case (W = 0.1290x10"®) the mini- 
mum film thickness is off to the sides in two areas close to the Hertzian circle. 

The variation of pressure and film thickness in the X-direction along a line close to 
the midplane of the conjunction is shown in figure 9 for three values of the dimensionless 
load parameter. The values of the dimensionless speed, material, and ellipticity pa- 
rameters were held fixed as described by equation (18) for all computations at various 
loads. In figure 9(a) as the dimensionless load is increased the inlet pressure becomes 
smaller. For the highest load case shown in figure 9(b), film thickness rises between 
the central region and the outlet restriction in the same manner as seen in figure 3(b) . 
Again this is attributed to the compressibility effects of the fluid. Also, at a load W of 
0. 5528x10"® the film thickness is slightly smaller than at a W of 1. 106x10“®. The 
reason is that the minimum film thickness is closer to the axial center of the contact at 
the lower load than at the higher load. As was pointed out in discussing figures 8(a) 
and (b), the location of the minimum -film-thickness region changes as the dimensionless 
load is changed . 


EFFECT OF MATERIAL PROPERTIES 

The effect of the dimensionless material parameter on minimum film thickness is 
not a simple matter. As can be seen from equations (3), (4), and (5), when either the 
material of the solids (as expressed in E’) or the material of the lubricant (as expressed 
in t] q and p iv ) is varied, not only does the material parameter G change, but so 
do the dimensionless speed U and load W parameters. Only the ellipticity parameter 
can be held fixed. For all the results presented in this section the ellipticity parameter 
is held fixed at 6. 

Table IV gives the four material -parameter results. The general form showing how 
the minimum film thickness is a function of the dimensionless material parameter is 
given as 


C = TG V 


(23) 


where 


C = 


H . 

min 

_ e “0. 68kjjj0. 68^-0. 073 


(24) 
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In equation (24) the exponents are rounded off to two significant figures so that any error 
could be absorbed in T given in equation (23) . By applying a least-square power fit to 
the four pairs of data, the values of T and V were found to be 


T = 3.6891 (25) 

V = 0.48669 « 0.49 (26) 

The coefficient of determination for these results was 0.9980, which is excellent. Sub- 
stituting equations (25) and (26) into equation (23) gives the values of H min shown in ta- 
ble IV. The percentage difference shown in table IV varies by less than 2 percent 
in all cases. Therefore, from equations (23) and (25) the effect of the dimensionless ma- 
terial parameter on the dimensionless film thickness can be written as 

<«> 


MINIMUM - FILM -THICKNESS FORMULA 


The proportionality expressions (11), (17), (22), and (27) established how the mini- 
mum film thickness varies with the ellipticity, speed, load, and material parameters, 
respectively. This enables a composite minimum -film -thickness formula for a fully 
flooded, isothermal, elastohydrodynamic point contact to be written as 


H . 
min 


3.63 u 0 - 68 G°- 49 W- 0 - 073 


(1-e 


-0. 68kj 


(28) 


In equation (28) the constant 3.63 is different from that in equation (25) to account for 

rounding off the material -parameter exponent. 

Table V gives the 34 different cases used in obtaining equation (28). In this table, 

H min corresponds to the minimum film thickness obtained from the EHL point-contact 

theory developed in reference 2, and H m ^ n is the minimum film thickness obtained from 

equation (28) . The percentage difference between these two values is expressed by D^, 

which is defined in equation (12) . In table V the values of are within ±5 percent. 

It is sometimes more convenient to express the side -leakage factor in equation (28) 

in terms of the radius of curvature ratio R /R instead of the ellipticity parameter k 

y x 

through the following relation: 

, R vO.64 

k = 1.03[ — ^ | (29) 

vv 
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where 


R y r Ay r By 

A = ^L + ^L 

R x r Ax r Bx 

Using equation (29) avoids the need to evaluate elliptic integrals of the first and second 
kinds in the determination of k. The minimum film thickness can thus be derived di- 
rectly from a knowledge of the radii of curvature of the contacting bodies (r^ x , r Bx , r^ y , 
and r By ) . 

It is interesting to compare the new, point -contact, minimum-film-thickness for- 
mula (eq. (28)) with the corresponding equation generated in the 1960’s (ref. 7) for line 
contacts 


(30) 

(31) 


R min,L 


= 2.65 u0.70q0.54w-0.13 


(32) 


The powers of U, G, and W in equations (28) and (32) are quite similar considering the 
different numerical procedures upon which they are based. It is also worth noting that 
the power of W in equation (28) is extremely close to the value of -0.074 proposed by 
Archard and Cowking in their study of point contacts (ref. 8). 


CENTRAL-FILM-THICKNESS FORMULA 

There is interest in knowning the central film thickness, in addition to the minimum 
'ilm thickness, in elastohydrodynamic contacts. The procedure used in obtaining the 
central film thickness was the same as that used in obtaining the minimum film thickness 
ind is not repeated here. The central -film-thickness formula obtained from the results 
s 


H c = 2.69 u 0<67 G 0 - 53 W" 0 ' 067 (1 - 0. 61 e“°- 73k ) (33) 

Comparing the central -film -thickness formula (eq. (33)) with the minimum -film- 
thickness formula (eq. (28)) reveals a slight difference. In equation (33) the load expo- 
nent is small but negative, as it was for the minimum -film -thickness formula. This is 
in contrast with the recent numerical study of Ranger, et al. (ref. 9) who found a small 
but positive exponent on the dimensionless load W in their formulation of a central film 
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thickness. 

Table VI gives the 34 different cases used to obtain equation (33) . In this table, H c 
corresponds to the central film thickness obtained from the EHL point-contact theory 
developed in reference 2, and H corresponds to the central film thickness obtained 
from equation (33) . The percentage difference between these two values is given by Dg 
and is written as 




(34) 


In table VI the values of Dg are within ±10 percent. 


CONCLUDING REMARKS 

By using the procedures outlined by the authors in an earlier publication the influ- 
ence of the ellipticity parameter and the dimensionless speed U, load W, and material 
G parameters on minimum film thickness has been investigated. The ellipticity param- 
eter was varied from 1 (a ball -on -plate configuration) to 8 (a configuration approaching 
a line contact) . The dimensionless speed parameter was varied over a range of nearly 
two orders of magnitude. The dimensionless load parameter was varied over a range of 
one order of magnitude. Situations equivalent to using solid materials of bronze, steel, 
and silicon nitride and lubricants of paraffinic and naphthenic mineral oils were con- 
sidered in investigating the role of the dimensionless material parameter. Thirty -four 
different cases were used to generate the minimum -film-thickness and central -film- 
thickness relations: 


H . = 3.63 u°* 68 G 0,49 W“ 0,073 (1 - e" 0- 68k ) 
min v ' 

H = 2.69 u 0 - 67 G 0,53 W 0 - 067 (1 - 0.61 e‘ 0 ' 73k ) 

V 


It was found that the ellipticity parameter k can be written as 


( 


k = 1.03 



where R /R is the radius of curvature ratio. 

y x 

Contour plots have been presented that indicate in detail the pressure distribution 
and the film thickness. In some solutions, pressure spikes were in evidence. The theo- 
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reticai solutions of. film thickness have all the essential features of previously reported 
experimental observations based upon optical Interferometry. 

The importance of the present report lies in the fact that it presents for the first 
time a complete theoretical film-thickness equation for elastohydrodynamic point con- 
tacts operating under fully flooded conditions. The exponents on the various dimension- 
less parameters governing minimum film thickness in such conjunctions are quite simi- 
lar to those developed earlier by Dowson for line contacts. The most dominant exponent 
occurs in association with the speed parameter, while the exponent on the load parameter 
is very small and negative . The material parameter also carries a significant exponent, 
although the range of this parameter in engineering situations is limited. A central- 
film -thickness formula exists for the contact geometry of a ball on a plate from which an 
estimate can be made of the minimum film thickness. However, the formula presented 
herein is valid for any contact geometry and proceeds directly to the evaluation of the 
minimum film thickness. 

Perhaps the most significant feature of the proposed minimum -film -thickness for- 
mula is that it can be applied to any contacting solids that present an elliptical Hertzian 
contact region. Many machine elements, particularly rolling -element bearings, poss- 
ess such geometry. And it is expected that the new minimum -film -thickness formula 
will find application in such fields. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 28, 1976, 

505-04. 
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TABLE I. - EFFECT OF ELLIPTICITY PARAMETER ON 


MINIMUM FILM THICKNESS 


Ellipticity 

parameter, 

k 

Minimum film thickness 

Difference between 

H min S min’ 

D P 

percent 

Obtained from 
EHL point- 
contact theory, 
H 

min 

Obtained from 
least- square 
fit, 

! no 

H 

"min 

i 

3. 367X10' 6 

3. 464X10" 6 

+2. 88 

1.25 

4. 105 

4.031 

-1.80 

1. 5 

4. 565 

4. 509 

-1. 22 

1.75 

4. 907 

4.913 

+. 11 

2 

5.255 

5. 252 

-.05 

2. 5 

5.755 

5.781 

+. 45 

3 

6.091 

6. 156 

+1.08 

4 

6.636 

6.613 

-.34 

6 

6. 969 

6.961 

-. 12 

8 

7.048 

7.050 

+. 02 


TABLE n. - EFFECT OF DIMENSIONLESS SPEED 
PARAMETER ON MINIMUM FILM THICKNESS 


Dimensionless 
speed param- 
eter, 

U 

Minimum film thickness 

Difference between 

no 

H min H min> 

D l> 

percent 

Obtained from 
EHL point - 
contact theory, 
H 

min 

Obtained from 
least-square 
fit, 

H 

min 

0. 08416X10’ 11 

3. 926xlO' 6 

3. 915X10" 6 

-0.275 

. 1683 

6. 156 

6.252 

+1. 564 

.2525 

8. 372 

8. 223 

-1.780 

.3367 

9. 995 

9.987 

-.078 

.4208 

11.61 

11.61 

-.004 

. 5892 

14. 39 

14. 57 

+1.280 

. 8416 

18. 34 

18. 54 

+1. 104 

1. 263 

24. 47 

24. 39 

-. 320 

1.683 

29.75 

29.61 

-.467 

2. 104 

34. 58 

34. 43 

-.432 

2. 525 

39.73 

38. 95 

-1.977 

2. 946 

43.47 

43. 22 

576 

3.367 

47.32 

47. 30 

-.042 

4.208 

54. 57 

54. 99 

+.765 

5.050 

61. 32 

62.20 

+1. 430 





TABLE m. - EFFECT OF DIMENSIONLESS LOAD PARAMETER 


ON MINIMUM FILM THICKNESS 


Dimensionless 
load param- 
eter, 

W 

Minimum film thickness 

Difference between 

rsj 

H min H min’ 

D l> 

percent 

Obtained from 
EHL point- 
contact theory, 

XT 

min 

Obtained from 
least- square 

H 

min 

0. 1106x10 " 6 

6. 969X10" 6 

6. 941X10' 6 

-0.41 

.2211 

6.492 

6. 599 

+1.65 

.3686 

6.317 

6.358 

+.64 

.5528 

6.268 

6. 172 

-1. 52 

.7371 

6. 156 

6.044 

-1. 81 

.9214 

6.085 

5. 947 

-2.27 

1. 106 

5. 811 

5. 868 

+.98 

1.290 

5.657 

5. 803 

+2. 58 


TABLE IV. - EFFECT OF SOLID MATERIAL AND LUBRICANT AS REPRESENTED IN DIMENSIONLESS 
MATERIAL PARAMETER ON MINIMUM FILM THICKNESS 


Solid ma- 

Lubricant 

Dimensionless 

Dimensionless 

Dimensionless 
load param- 
eter, 

W 

Minimum film thickness 

terial 


material pa- 
rameter, 

G 

speed param- 
eter, 

U 

Obtained from 
EHL point - 
contact theory, 

ii 

min 

Obtained from 
least- square 
fit, 

5 min 

Bronze 

Paraffinic 

2310 

0. 3296X10" 11 

0. 7216xl0“ 6 

6.931X10" 6 

6. 873X10 -6 

Bronze 

Naphthenic 

3591 

.9422 

.7216 

17. 19 

17. 404 

Steel 

Paraffinic 

4522 

. 1683 

.3686 

6.317 

6.336 

Silicon 

nitride 

Paraffinic 

6785 

. 1122 

.2456 

6.080 

6.038 


Difference between 
H min 311(1 S min’ 


percent 


-0. 84 
+1.25 
+.31 
-.70 
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TABLE V. - DATA SHOWING EFFECT OF ELLIPTICITY, LOAD, SPEED, AND MATERIAL ON MINIMUM FILM THICKNESS 


Case 

Eliipticity 

Dimensionless 

Dimensionless 

Dimensionless 

Minimum film thickness 

Difference between 

Results 


parameter, 

load param- 

speed param- 

material pa- 



H min and ”min’ 




k 

eter. 

eter. 

meter, 

Obtained from 

Obtained from 

D P 






W 


U 

G 

EHL point- 

least-square 

percent 










contact theory, 

fit, 











H min 

H min 



1 

1 

0. 1106X10 -6 

0. 1683X10" 11 

4522 

3. 367X10" 6 

3. 514X10" 6 

+4.37 



2 

1.25 







4. 105 

4.078 

-.66 



3 

1.5 







4. 565 

4. 554 

-. 24 



4 

1.75 







4. 907 

4.955 

+.98 



5 









5.255 

5.294 

+.74 


^Eliipticity 

6 


1.5 







5.755 

5. 821 

+1. 15 



7 

3 







6. 091 

6. 196 

+1.72 



8 

4 







6. 636 

6.652 

+. 24 



9 

« 








6.969 

7. 001 

+. 46 



10 

8 







7. 048 

7.091 

+.61 



11 

6 

.2211 





6.492 

6.656 

+2.53 



12 



. 3686 





6. 317 

6. 412 

+1. 50 



13 



. 5528 





6. 268 

6.225 

-.69 


Load plus 

14 



.7371 





6. 156 

6.095 

-.99 


case 9 

15 



.9214 





6.085 

5.997 

-1.45 



16 



1. 106 





5. 811 

5.918 

+1. 84 



17 



1.290 





5.657 

5.851 

+3.43 

- 


18 



.7371 

.08 

416 



3.926 

3. 805 

-3.08 



19 





.2525 



8. 372 

8.032 

-4.06 



! 20 





.3367 



9.995 

9.769 

-2.26 



: 21 





.4208 



11.61 

11. 37 

-2.07 



22 





.5892 



14. 39 

14. 29 

-.69 



j 23 





. 8416 



18.34 

18.21 

-.71 


Speed plus 

24 





1.263 



24. 47 

24.00 

-1.92 


^ case 14 

25 





1.683 



29.75 

29. 18 

-1. 92 



26 





2. 104 



34. 58 

33.96 

-1.79 



27 





2. 525 



39. 73 

38.44 

-3.25 



28 





2.946 



43. 47 

42.69 

-1.79 



29 





3. 367 



47. 32 

46.76 

-1. 18 



30 





4. 208 



54. 57 

54.41 

-.29 



31 





5.050 



61. 32 

61. 59 

+.44 



32 



.7216 

. 3296 

2310 

6.931 

6.938 

+. 10 

''j Materials 

33 



.7216 

. 9422 

3491 

17. 19 

17. 59 

+2. 33 

> plus 

34 



.2456 

. 1122 

6785 

6. 080 

6. 116 

+. 59 

J case 9 



TABLE VI. - DATA SHOWING EFFECT OF ELLLPTICITY, LOAD, SPEED, AND MATERIAL ON CENTRAL FILM THICKNESS 


Case 

Ellipticity 

Dimensionless 

Dimensionless 

Dimensionless 

| Central film thickness 


parameter, 

k 

load param- 
eter, 

W 

speed param- 
eter, 

U 

material pa- 
rameter, 

G 

Obtained from 
EHL point- 

Obtained from 
least-square 








contact theory, 

fit, 








H c 

«c 

1 

1 

0. 1106X10’ 6 

0.1683x10' 11 

4522 

6. 860xl0~ 6 

6. 215X10" 6 

2 

1.25 







6.964 

6.647 

3 

1. 5 







7.001 

7.006 

4 

1.75 







7.015 

7. 306 

5 

2 







7.402 

7. 556 

6 

2.5 







7.653 

7. 937 

7 

3 







7. 845 

8. 202 

8 

4 







8.292 

8. 513 

9 

6 







8.657 

8. 736 

10 

8 







8.672 

8. 787 

11 

6 

.2211 





7.796 

8.339 

12 



. 3686 





7. 505 

8. 059 

13 



. 5528 





7. 309 

7. 843 

14 



.7371 





7. 517 

7.693 

15 



. 9214 





7.611 

7. 578 

16 



1.106 





7.416 

7. 487 

17 



1.290 





6.762 

7.410 

18 



. 7371 

.08416 



4.917 

4. 836 

19 





.2525 



9. 999 

10. 10 

20 





. 3367 



11.40 

12. 24 

21 





. 4208 



13.07 

14.21 

22 





. 5892 



17. 13 

17. 81 

23 





. 8416 



21.35 

22.61 

24 





1.263 



29.62 

29.68 

25 





1.683 



35. 50 

35. 98 

26 





2. 104 



41.05 

41.79 

27 





2. 525 



46.64 

47.22 

28 





2. 946 



51.08 

52. 36 

29 





3. 367 



55. 56 

57. 26 

30 





4. 208 



63.81 

66.49 

31 





5.050 



71.25 

75. 13 

32 



.7216 

. 3296 

2310 

8. 422 

8. 466 

33 



.7216 

. 9422 

3591 

20.70 

21.62 

34 

T ! 

.2456 

. 1122 

6785 

7. 825 

7. 825 


Difference between 


2 ’ 


percent 


Results 


-9.40 
-4. 55 
+. 07 
+4. 15 
+2.08 
+3.71 
+4. 55 
+2.67 
+.91 
+1.33 
+6.97 
+7.38 
+7. 31 
+2.34 
-.43 
+. 96 
+9. 58 
-1. 65 
+ 1.01 
+7.37 
+8.72 
+3. 97 
+5. 90 
+.20 
+1. 35 
+1. 80 
+1.24 
+2. 51 
+3.06 
+4.20 
+5. 45 
: . 52 
+4.44 
0 


/-Ellipticity 


ULoad plus 
I case 9 


J 


I Speed plus 
f case 14 


"'j Materials 
i plus 
J case 9 











Dimensionless pressure, P = p/E' 




1.25 



(b) Dimensionless film thickness. 
Figure 3. - Concluded. 
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Dimensionless 
film thickness, 
H*h/R x 

A 61.4xl0' 6 
B 62.0 
C 63.0 
D 65.0 
E 68.0 
F 72.0 
G 78.0 
H 86.0 
H 



(b) U = 0.5050X10' 10 . 
Figure 5. - Concluded. 
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Dimensionless pressure, P = p/E' 



(a) Dimensionless film thickness. 

Figure 6. - Variation of dimensionless pressure and film thickness on X-axis for 
three values of dimensionless speed parameter. The value of Y is held fixed 
near axial center of contact. 
















Dimensionless 
film thickness, 
H = h/R x 

A 7. 0x10” 6 

B 7.2 

C 7.4 

D 7.7 

E 8.1 

F 8.7 

G 9.5 

H 10.7 



(a) W ■ 0.1106xl0’ 6 . 

Figure 8. - Contour plots of dimensionless film thickness for dimensionless load 
parameters W of 0.1106x10’° and 0.1290xl0" 5 , respectively. The dimensionless 
parameters k, W, and G are held constant as defined in equation (18). 




ia» Dimensionless pressure. 


Figure 9. - Variation of dimensionless pressure and film thickness on X-axis for 
three values of dimensionless load parameter. The value of Y is held fixed near 
axial center of contact. 










